arXiv: 1503.07814vl [math-ph] 26 Mar 2015 


Perturbative Construction of Models of Algebraic 
Quantum Field Theory 


Klaus Fredenhagenl 1 ), Katarzyna Rejzneb 2 ) 


C) II Institute for Theoretical Physik, University of Hamburg, 
Department of Mathematics, University of York 
klaus.fredenhagenSdesy.de, kasia.rej znerOyork.ac.uk 

March 27, 2015 


Abstract 

We review the construction of models of algebraic quantum field theory by renor¬ 
malized perturbation theory. 


1 Introduction 


The axiomatic framework of AQFT allows for a qualitative description of a large 
class of phenomena occurring in particle physics and some parts of solid state physics. 
It does not, however, yield quantitative predictions, and there is a widespread im¬ 
pression that one has to abandon the formalism of AQFT if one wants to make real 
contact with experiments. Actually, up to now no single model of an interacting 
AQFT in 4d Minkowski space has been constructed. 

But what are the alternatives? Standard textbooks on QFT either start from 
canonical quantization of free field theory on Fock space and try to construct the 
interacting theory in the interaction picture, or they use the path integral formalism. 
The canonical approach ends up in the Gell-Mann Low formula for the vacuum 
expectation values of time ordered products of fields, 


uo{Tip{x\)... <p(x n )) 


(U, Ty 0 (a:i)... ip 0 (x n )eK 7 £r(U d Vn) 
(H,TeA f 


( 1 ) 


where tpo is the free field treated as an operator valued distribution on the Fock 
space, Lj is the interaction density treated as a Wick polynomial of tpo and is 
the vacuum vector of the free theory. The time ordering symbol T means that the 
products have to be performed after ordering of the factors according to their time 
arguments. 
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The path integral approach reinterprets the Gell-Mann Low formula as an inte¬ 
gral over all classical field configurations <t> 

u 0 (T(p(xi)...ip(x n )) = Zr'j (t>(xi).. .(f>{x n )e^ ^ c{x)dix Dcj) (2) 

where now L is the full classical Lagrangian, Z is a normalization factor, and Dcj) is 
thought of as the Lebesgue integral over field space. 

Both versions are only heuristic, and it required the hard and ingenious work of 
several generations of physicists to turn these formal expressions into unambiguous 
computations. The state of the art is that one can create a formal power series in h 
where every term is well defined, up to some remaining infrared problems originating 
from the integral over Minkowski space in the exponent. The great success of QFT 
relies on the fact that already the first few terms of this series yield a good and often 
even excellent agreement with experimental data. 

The path integral approach has the advantage that it is formally similar to 
probability theory. Actually, by passing to imaginary time (Wick rotation), one 
can interpret the vacuum expectation values of time ordered products of fields as 
correlation functions of a probability distribution (euclidean QFT). In particular, 
counter-intuitive properties of quantum physics as e.g. entanglement do not occur. 
Moreover, the momentum space integrals in the evaluation of Feynman diagrams 
have better convergence properties. Finally, due to the Osterwalder-Schrader theo¬ 
rem, a Wick rotation back to real time is possible under very general conditions. 

The disadvantage of the path integral approach is that the nonconrmutative 
product of operators, which is crucial for the structure of quantum physics, appears 
only indirectly in terms of different boundary values of analytic functions. In the 
canonical approach, the operator product is given from the beginning, but there 
the definition of the time ordered product is problematic. First of all, it is not well 
defined as a product of operators, since, by the existence of a deterministic time 
evolution, fields at a given time can be expressed in terms of fields at an earlier time, 
and thus the time ordering prescription is ambiguous. One may instead define time 
ordered products TA(t.\)... A(t n ) of an operator valued function of time t. A(t) 
as a symmetric operator valued function of n time variables such that 


TA(ti)... A(t n ) = A(ti)... A(t n ) if h> ■■■ >t n . (3) 

This, however, does not work since the quantum fields are distributions, and the 
time ordering prescription would amount to multiply them with a discontinuous 
function. 

But there is a way out, as first observed by Stiickelberg, further elaborated by 
Bogoliubov and collaborators and finally worked out by Epstein and Glaser (causal 
perturbation theory). Namely, one may define the time ordered product of n fields 
as an operator valued distribution which is already known for non-coinciding points. 
Due to the UV divergences of QFT, the extension to coinciding points is ambiguous, 
but the crucial observation is that this ambiguity is the same ambiguity which occurs 
in the removal of infinities in approaches where the theory is regularized by the 
introduction of a momentum cutoff, and where the theory without cutoff has to be 
fixed by renormalization conditions. 

Originally, the insertion of a test function g into the interaction Lagrangian, 
instead of integrating it over all spacetime, was considered to be an intermediate 
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step, and in the last step one aimed at the limit where g tends to 1 (adiabatic 
limit). In this limit one then finds vacuum expectation values of operator products of 
time ordered products of interacting fields, and using the Wightman reconstruction 
theorem, one obtains interacting fields as operator valued distributions on some 
“Hilbert space”, of course only in the sense of formal power series. But as first 
observed in [32] and rediscovered in [5j the algebra of observables associated to 
some bounded region can be already constructed if one chooses a test function g 
which is equal to 1 on some slightly larger region. Actually, the full Haag-Kastler 
net of the interacting theory can be obtained in this way. Thus causal perturbation 
theory provides a direct way for a construction of the algebra of observables. Hence 
by replacing the condition that the local algebras have to be unital C*-algebras by 
the condition that they are isomorphic to unital *-algebras of formal power series of 
operators on a dense invariant subspace of some Hilbert space, one obtains a huge 
class of models, in particular the models used in elementary particle physics. 

On this level, structural properties of the local net can be analyzed, but the 
powerful structural results on C*- and von Neumann algebras are not available. 
Nevertheless, one can derive interesting results, as e.g. the validity of the time-slice 
axiom nu, the existence of operator expansions [25] and an algebraic version of the 
Callan-Symanzik equation [5]. 

In order to reach numerical predictions, one needs in the next step a construction 
of states. States are here defined as linear maps from the algebra to the formal power 
series over C, and the positivity condition on states now means that the expectation 
value of A* A is the absolute square of another power series. The construction of 
states can be done via the adiabatic limit as described above; this is the way the 
vacuum state is constructed in m • As observed by Steinmann BE this method does 
not work for the construction of KMS states. The reason is that the analog of the 
Gell-Mann Low formula does not hold at nonzero temperature, due to the different 
asymptotic time behavior of free and interacting systems at nonzero temperature. 
But here a structural result helps: namely, the time-slice axiom allows to treat only 
the theory within a short time interval, and the asymptotic behavior in time does 
not matter for the existence of states. What matters is the decay of correlations in 
spacelike directions which is exponentially fast for massive theories. 

Up to now we considered the so-called on shell formalism, due to the fact that 
we constructed the operators on Fock space, thereby imposing the validity of the 
Klein Gordon equation for the free field. It turned out, however, to be more useful to 
replace Fock space operators by functionals of classical field configurations which are 
not restricted to those which satisfy the field equation. On the space of functionals 
one can then introduce several operations: the pointwise (classical) product, the 
involution by complex conjugation, the Peierls bracket (as a covariant version of 
a Poisson bracket on the space of functionals), the non-commutative, associative 
★-product in the sense of deformation quantization and the time ordered product. 
It is the latter which is relevant for inducing the interaction and which requires 
renormalization. The other operations can be directly defined. This is trivial for 
the pointwise product and for the involution. The Peierls bracket is obtained by 
considering the linearized Euler-Lagrange operator, which e.g. for the </? 4 -theory 
looks like 

E'(v) = n + m 2 + ^cp 2 (4) 

where the last term acts as a multiplication operator. We consider only theories 
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where the linearized Euler-Lagrange operator is normally hyperbolic and hence has 
unique retarded and advanced propagators A R t A (ip). The Peierls bracket is then 
defined by 

= ( 5 ) 

where A = A R — A" 4 and j- is the functional derivative (defined as the directional 
derivative). In free theories, E' and then also the propagators do not depend on cp. 
One then can define the ^-product (in the sense of formal power series in K) by 

(F*G)(ip) = e*\& A * r /F(<p)G(<f/) ■ (6) 

The time-ordered product is defined by a similar formula 

(F -T G) = e h (^’ A ° fc)F(<p)G(tp') (7) 

<p f =(p 

with the Dirac propagator A D = |( A R + A A ). 

Note that there is a crucial difference between the time ordered product and the 
other products. Namely, the ideal generated by the field equation with respect to 
the pointwise product is also an ideal with respect to the Poisson bracket and the 
★-product, but not with respect to the time ordered product. This is actually a 
necessary condition which allows to use the time ordered product for introducing an 
interaction. Let V be the interaction. We then define the interacting observables by 

R v (F) = (e%y- 1 *(eV-rF) . ( 8 ) 

Here eg- means the exponential series where powers are computed via the time 
ordered product. For an evaluation functional d’ a; (v ? ) = ( Pi x )i the corresponding 
interacting held x —> Ryi^x) satisfies the equation 

E'R V (<S> X ) = E'§ X + R V ( 9 ) 

which may be interpreted as the held equation with interaction V, when evaluated 
on some <p which satishes the free held equation. 

The rough description of the formalism has to be made precise in the following 
sense: One has to specify the functionals which are allowed, and one has to check 
whether this class contains the relevant ones. As it stands we need functionals whose 
functional derivatives are test functions in order that all operations are well defined. 
But we will see that by changing the products to equivalent ones which corresponds 
to Wick ordering in the Fock space framework one can extend the Peierls bracket 
and the ★-product to a rather large class of functionals, which contains in particular 
the local functionals that appear as terms in the Lagrangian, and is stable under 
these operations. The definition of time ordered products is more involved and there 
are different possibilities, corresponding to the choice of renormalization conditions. 

The plan of the paper is as follows: we will first outline the functional analytic 
tools which are needed for the operations. We then define the Peierls bracket and 
the ★-product on a class of functionals called microcausal. Thereafter we come to 
the problem to define the time ordered products. Here we first develop the general 
formalism and show that it leads to a construction of local nets. We review some 
structural properties of these nets, in particular their behavior under renormalization 
group transformations. Finally, we outline a possible construction of states. 
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2 Functional derivatives, wave front sets, and all 
that 


Our approach to quantum field theory is based, in this respect similar to the path 
integral approach, on functionals of classical field configurations. But there, at least 
in its euclidean version, the measure theoretic aspects of the space of field configu¬ 
rations are of central importance; in our case, due to the frequent use of functional 
derivatives, the properties of the space of field configurations as a differential mani¬ 
fold are crucial. 

For definiteness we concentrate on the case of a scalar field and fix an oriented, 
time-oriented globally hyperbolic spacetime M. There we consider the space of real 
valued smooth functions as the space of field configurations, 

£=e°°(Af,R). (10) 

We model it as a differentiable manifold over the space of compactly supported 
smooth functions 

D = e“(M,M) ( 11 ) 

where charts are defined as maps 

tp + "D “D, tp + tp >—¥ tp , with tp E £ . (12) 

Clearly, £ has the structure of an affine manifold. A similar affine structure can be 
introduced also for other fields, as e.g. gauge theories or gravity. 

We will model observables as functionals on £ and we allow only functionals 
which depend on the field configuration inside some compact region. This includes 
in particular the polynomial functionals 

n . 

n<p) = E / <p(xi)...ip(xk)fk(xi,---,Xk) ( 13 ) 

fc=l 

with symmetrical distributional densities fk with compact support. More generally, 
we consider functionals, for which all functional derivatives F (") exist and are con¬ 
tinuous. We recall after [24] (see |37] for a review) that a functional derivative of a 
functional is defined as 

{F {1) (p),p) :=-^F(p + A<£)| a=0 (14) 

and a functional is differentiable if the derivative exists for all p £ £. It is continu¬ 
ously differentiable if the map £ x T> —¥ C, (tp, tp) i-> {F^\tp), tp) is continuous. If £ 
is taken with its natural Frechet topology, this implies that F W [tp) is a compactly 
supported distributional density. Higher derivatives are obtained by iterating this 
definition, i.e. 

F {n) [tp)[0i,...,tp n ) ■■= ^x p( - n ~ l) [tp + Xipn)^,... ,tp n _i)\ x=Q , (15) 

and we find that the F^[p)’s are symmetric compactly supported distributional 
densities with support contained in K n , for some compact set K C M. There re¬ 
mains, however, the problem that the propagators have singularities, and therefore 
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the contractions with the distributional densities occurring as functional deriva¬ 
tives are not always well defined. The restriction to functionals whose functional 
derivatives are smooth densities, on the other side, would exclude almost all local 
functionals, i.e. functionals of the form 

F (<P) = j fUx(<P)) , (16) 

where jx{'~p) = (x, <f{x), dip(x ),...) is the jet prolongation of <p and / is a density¬ 
valued function on the jet bundle. For these functionals, the derivatives are sup¬ 
ported on the thin diagonal 

D n = {(*i) •• ■ ,x„) G M n ,x i = • • • = x n } (17) 

and thus smooth for n > 1 only when they vanish. 

The singularities of distributions can be analyzed using the concept of the wave 
front set. On Minkowski spacetime, this concept arises in the study of the decay 
properties of the Fourier transform of the given distribution multiplied by some test 
function. A pair of a spacetime point x and a nonzero momentum k is an element of 
the wave front set of a given distribution t if, for any test function / with /( x) ^ 0 
and some open cone around k, the Fourier transform of ft does not decay fast (i.e. 
faster than any power) inside the cone. The notion of the WF set can be generalized 
to an arbitrary smooth manifold M, and it is defined as a subset of T*M. 

As the first example we will consider the Dirac S distribution. Since (fS, e lkm ) = 
/(0), where e lfe *(x) = e lkx , if follows that for any choice of / with f(x) ^ 0 the 
Fourier transform of ft does not decay fast in any direction and hence the wave 
front set of S is 

WF(<5) = {(0, k), k 0} . (18) 

Another important example is the distribution / >->• lim e to f i+^dx. Note that its 
Fourier transform is 


lim 

eto 


/ 


f (x) A l, T T 

e lkx dx = -* 

x + ie 



f{k')dk' . 


and / A °° f{k')dk' decays strongly as k —> oo, while for k —> —oo we obtain 


lim 

k — y —oo 



f(k')dk' 


2tt/(0) . 


(19) 


( 20 ) 


We can now conclude that 

WF(lim(x + ie) -1 ) = {(0, fc), k < 0} . (21) 

e4,0 

For more information on wave front sets see m or chapter 4 of [ 3 ] . Using WF 
sets we can formulate a sufficient condition for a pointwise product of distributions 
to be well defined. Let t and s be distributions on M. The Whitney sum (i.e 
pointwise sum) of their wave front sets is defined by 

WF(f) + WF(s) = {(x, k + k') |(x, k) G WF(f), (x, k') G WF(s)} (22) 
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If this set does not intersect the zero section of T*M, then we can define the pointwise 
product ts as 

(ts, fg) = J tf(k)sg(—k)dk , (23) 

where /, g £ T> are chosen with sufficiently small support. To see that the integral 
above converges, note that if k ^ 0 , then either tf is fast-decaying in a conical 
neighborhood around k or sg is fast-decaying in a conical neighborhood around —k, 
while the other factor is polynomially bounded. 

Beside the criterion for multiplying distributions, WF sets provide also a charac¬ 
terization of the propagation of singularities. Let P be a partial differential operator 
and ap its principal symbol. We can interpret ap as a function on the cotangent 
bundle T*M : which carries a structure of a symplectic manifold. With the use of the 
canonical symplectic form, 1-forms on T*M can be canonically identified with vector 
fields. Let Xp be the vector field (called the Hamiltonian vector field) corresponding 
to the 1-form dap. In coordinates it is given by 

^ dap d dap d 

P i^i 9xj dxj dkj 


Let (xj(t),kj(t)) be a curve that fulfills the system of equations (Hamilton’s equa¬ 
tions) : 

dxj dap 
dt dkj 1 
dkj dap 

dt dxj 

We call a solution (xj(t), kj(t )) of the above equations an integral curve of Xp and 
the bicharacteristic flow is defined as the set of all such solutions. Along this flow 
f= Xp(ap) = 0, so ap is conserved under the bicharacteristic flow. We are now 
ready to state the theorem on the propagation of singularities: the wave front set 
of a solution u of the equation Pu = f with / smooth is a union of orbits of the 
Hamiltonian flow Xp on the characteristics char P = {(a;, k ) £ T* M\a{P)(x, k ) = 0} 
of P. 

For hyperbolic differential operators on globally hyperbolic spacetimes (for ex¬ 
ample M), the set of characteristics is the light cone, and the principal symbol is 
the metric on the cotangent bundle. For such operators the wave front set of solu¬ 
tions is therefore a union of null geodesics 7 together with their cotangent vectors 

k = g{ 7 , ■)• 


3 The Peierls bracket and the ^-product 

As outlined in the Introduction, we start our construction of a pAQFT model from 
the classical theory. To this end, we equip the space of functionals on the config¬ 
uration space with a Poisson structure provided with the so called Peierls bracket. 
This bracket, introduced in [3B], is the off-shell extension of the canonical bracket of 
classical mechanics, which is defined only on the space of solutions to the equa¬ 
tions of motion. To see how this works, we will start in a setting which resembles 
closely classical mechanics and then show the relation with the Peierls method on a 
concrete example. 
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3.1 Canonical formalism and the approach of Peierls 

Let us start with the free scalar field with the field equation 

P<P — 0 , ( 24 ) 

where P = □ + to 2 is the Klein-Gordon operator. For this equation the retarded 
and advanced Green’s functions exist. We also know that for every / £ £ whose 
support is past and future compact, A/ is a solution to ll2Hl . Conversely, every 
smooth solution of the Klein Gordon equation is of the form A/ for some / £ £(M) 
with future and past compact support. 

Without loss of generality, the spacetime can be assumed to be of the form 
M = R x E with Cauchy surfaces {t} x E, t £ R. The space of Cauchy data 
(ifi(t, x ), tp(t, x)) on the surface {t} x E is 

e = {(</>,*/>) £ £(E) x £(E)}, 

where £(E) = C°°(E,R). This space is isomorphic to £g, the of smooth solutions to 

(EH)- 

As in classical mechanics, equations of motion can be derived from the least ac¬ 
tion principle. Elements of C play the role of generalized coordinates and generalized 
velocities, while a smooth trajectory t i—</>(f), t £ R is a function which assigns to 
an instant of time t a function <j>(t) £ £(E) such that trajectories </> are in one to one 
correspondence with field configurations ip : (t,x) —> <j>(t)(x), i.e. elements of £. 

The Lagrangian L associates to every compact region KcEa functional Lx 
on C, typically given in terms of a Lagrangian density L, 

L K {(/),ip)= [ £(<t>(x), V4>{x),ip(x))d<r(x) , 

Jk 

and the action is, for every compact K C E and every finite time interval I, a 
function on the space of trajectories defined by 

Sixk{4>) = J L K {<j>{t),4>{t))dt = J (/ £(<p(t,x), V x (p(t,x),ip(t,x))da t (x) S j dtt . 

.. _ ( 25 ) 

Solutions are configurations for which, for all compact K and /, Sj x k is stationary 
under variations 5(j> with support in the interior of I x K. If e.g. L is the Lagrangian 
density of the free scalar field, then the least action principle yields (1241) as the 
equation of motion. 

Now let F, G be two functions on the space of trajectories which depend only on 
the restriction of the trajectory to [ti, t^\ x K for some compact K C E and t\ < t?. 
Let £s be the space of solutions for an action S, and let t\g : £s —> £s+ag be 
the map which associates to a solution for S a solution for S + A G such that both 
solutions coincide for t <t\ {r\o is called the retarded Mqller map). Following the 
idea of Peierls, we consider the change of F under the change of the action and set, 
for a solution </> £ £ 5 , 

D gF{ 4>) = -j^\\=oF(r\ G (<t>)) ■ 

Similarly, we introduce the advanced Mqller map ciaf : £s —>■ &s+\f where the 
solutions coincide for t > t- 2 , and set 

GfG(</>) = -tv|a=oG(oaf(</>)) . 


8 


The Peierls bracket of G and F is now defined by 


{G,F} Pei = B G F-a F G. (26) 

The advantage of the Peierls bracket is the fact that it is defined covariantly, directly 
in the Lagrangian formalism. As it stands, the Peierls bracket of two functionals is 
only defined on solutions, and one has to prove that it depends only on the restriction 
of the functionals to the space of solutions. In order to show that it satisfies the 
Jacobi identity, one has to extend it to a neighborhood of the space of solutions. It is, 
however, possible and also convenient to extend it to a Poisson bracket on functions 
of arbitrary configurations £ (not only of solutions £g). We will now derive another 
formula for the Peierls bracket (formula (EJ) from the Introduction), which makes 
use of retarded and advanced Green’s functions for normally hyperbolic operators. 
Next we will show that, restricted to the solution space, flSJ) is equivalent to the 
canonical bracket. 

3.2 The generalized Lagrangian formalism 

Before we continue, there is one small modification to the classical Lagrangian for¬ 
malism, which we have to perform in order to make the quantization simpler. In 
formula ([25)l . we have smeared the Lagrangian density L(x,t) with a characteristic 
function of a certain compact region. Such sharp cut-offs would introduce additional 
divergences in the quantum theory, which we wish to avoid. Therefore, we replace 
the characteristic function by a smooth function that is equal to 1 on a sufficiently 
large region. Actually, it is convenient to consider all possible cutoffs and define the 
generalized Lagrangian as a map L from T) to the space 3J oc of local functionals on 
£. We require that 


Mf + 9 + h) — L{f + g) — L{g) + L(g + h ), 
for /, < 7 , h £ D and supp / (~l supp h = 0 . We also want 

supp(L(/)) C supp(/), 

where the support of a smooth functional F £ C°°(£,C) is defined as 

supp F = {x £ M|V neighborhoods U of x 3 ip, ip £ £, supp ip C U , (27) 

such that F(tp + ip) ^ F(</j)} • 

The action is an equivalence class of Lagrangians, where L\ ~ L 2 if 

supp(Li - L 2 ){f) C supp df . (28) 

The Euler-Lagrange derivative is a map S' : £ —> D' defined as 

<S'(^),h) = (L(/)«( V ),/i) , (29) 

with / = 1 on supp/i. Note that S' £ T(T*£). The field equation is now the 
condition that 

S\<p) = 0, (30) 
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which coincides with the condition obtained from the variation of (ill). We model 
observables as multilocal functionals on £ (i.e. products of local functionals). The 
maps F, G considered in the previous section are examples of such functionals. 
The space of multilocal functionals on the space of solutions to m is given by the 
quotient T/To, where To denotes the space of multilocal functionals that vanish on 
£s- 

The second variational derivative of the action is defined by 
(S"(p), hi <g) h 2 ) = ^ L {2 \f)(p ), hi <g> h^J , 


where / = 1 on supp hi and supp h 2 . S" defined in such a way is a symmetric 
two tensor on the affine manifold £ (equipped with the smooth structure induced 
by tw) and for each p it induces an operator from D to D'. Moreover, since L(f) is 
local, the second derivative has support on the diagonal, so S"(p) can be evaluated 
on smooth functions hi, h 2 , where only one of them is required to be compactly 
supported. This way we obtain an operator (the so called linearized Euler-Lagrange 
operator) E'[S](p) : £ —>■ T>'. 

We want to show now that the original formula of Peierls (1261) is equivalent to 
(|5j). if E'[S](p) is a normally hyperbolic operator. Let G G Ti oc be a local functional. 
We are interested in the flow (4 >a) on £ which deforms solutions of the original field 
equation S'(ip) = 0 to those of the perturbed equation S'(p) + A G^\p) = 0. Let 
$o(y>) = P and 


(s' M ($ x (p))+GW(3> x (p)j) 


= 0 . 

A—0 


(31) 


The vector field p t->- X(p) = ^$a( < / : ’)|a=o satisfies the equation 


(E'[S](p),X(p)) + G^(p) = 0. 


(32) 


RIA 

Let A s (ip) be the retarded/advanced Green’s function of the normally hyperbolic 
operator E'[S]((p ) and let A s(p) = A^(<^) — Ag(<p) be the causal propagator. We 
obtain now two distinguished solutions to the equation (132D . 


X R/A (P ) = (Af /A M, G& M) ■ (33) 


Note that X R (tp) = (Dg$) (tp), where $ is the evaluation functional <&x(p) = p(x). 
The difference X = X R — X A defines a vector field X G T(T£(M)) and it follows 
that 

{G,F} Pei (p) = D G F(p) - a F G(p) = (F^(p),A RIA (p)G^(p)) . 


Now we prove the equivalence between (0 and the canonical bracket. We fix a 
Cauchy surface {f} x E. Note that, given Cauchy data (<j>,ij)) G C, we can write the 
unique solution p corresponding to these Cauchy data as 

P(x) = P((/>,i(>)(x) = j ^A s (x;t,y)ip(y) - ^A s (x\t,y)(j)(y)^ da t (y). (34) 

Canonical momenta are obtained as distributional densities by 

(Tr((j),il)),h) = -^-|a=o L K ((/>,ip + A h), h G T)(K) . 
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We assume that for the Lagrangians of interest ir is always smooth. The phase space 
is then 

?=£(£) x £ d (E), (35) 

where £<*(£) is the space of smooth densities. The tangent space of 7 at 

some point (</>, if>) consists of the compactly supported elements (/, g ) G 7. The 
phase space has the canonical symplectic form 

CT(0»((/l,/ 2 ),(ffl,ff2)) = / (/1S2 -/2Sl)- 

JY 

Note that £(E) x £^(£) C £(£) x CD'(E) = T*(£(£)), so (CP, cr) is indeed the analog 
of the phase space in classical mechanics. 

For simplicity we consider an action S induced by a Lagrangian L which depends 
on cj) only through the kinetic term \(j) 2 , hence n(y) = </ )(y)da t (y ). Let a : (</>, n) >—>■ 
(4>, <j>) and /3 = /3oa:CP^-£s- We can now prove the equivalence of the canonical 
and the Peierls bracket. Let F, G € 7. Using (l34l) we obtain 

{Fo/3,Go/3} can = ((f§ ° s^j) (if ° M®)) ~ (if ° JF(x )) (if 0 

(e,£Wo^®G (1) o^ , 

where 0 is given by 

0(z',z) = J (A s (z';t,x)A s {z;t,x) - A s (z;t,x)A s (z'-,t,x) S ) da(x). 

From general properties of the causal propagator A g (the generalization of (1341) to 
distributional Cauchy data) it follows that the convolution 0 above is equal to As- 
Hence, on the solution space £ 5 , 

{foftGo ^} can = {F, G}pei o j3. 

3.3 Example: the Poisson bracket of the interaction 

In this section, following |21] . we give another argument for the equivalence of the 
Peierls and the canonical bracket on the example of the ip A interaction. Consider 
the generalized Lagrangian 

L (f)(v) = J ^v^V'V - fdy , 

M 

where d/j,(x) is the invariant measure fi on M, induced by the metric. Then S'((f) = 
-((□ + m 2 )<p + l)^ 3 ) and E'[S](ip) is the linear operator 

- (d + m 2 + ^(p 2 ^ (36) 

(the last term is to be understood as a multiplication operator). The Peierls bracket 
is 

{^,^}pei = A s ($)(a:,y), (37) 
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where $ 3 ;, are evaluation functionals on £ and x 1 —> As(p)(x, y) is a solution (at 
p) of the linearized equation of motion with the initial conditions 

AsM(y°, x; y°, y) = 0 , A s {p)(y°, x; y) = <5(x, y) . (38) 

This coincides with the Poisson bracket in the canonical formalism. Namely, let 
p £ £5, then 

0= j(n + m 2 )$, + ^,$ i/ J = (n + m 2 + ^ 2 ^j {$x,$J can • (39) 

In the first step we used the fact that p is a solution of the equations of motion and 
in the second step we used the fact that the canonical bracket is a derivation in both 
arguments. We can see from the equation above that the canonical Poisson bracket 
satisfies the linearized field equation with the same initial conditions as the Peierls 
bracket. The uniqueness of solutions to these linearized equations implies that in 
fact ^y}can = y }p e i, on £g. This clearly extends to general functionals 

since {P, G} Pei = ({$*,, $ y }pei, P (1) < 8 > G (1) ), and similarly for {., ,} can . 

3.4 Geometrical structures in classical theory 

The Peierls bracket, from now on denoted by {., .}s, introduces a symplectic struc¬ 
ture on the space T/To of on-shell multilocal functionals on £g. We can find a 
nice geometrical interpretation for this space using some basic notions of symplectic 
geometry. Let us assume that S is quadratic, so the equations of motion are of the 
form S'(p) = Pip = 0 for some normally hyperbolic differential operator P. The 
space of solutions £g is a vector space and hence an infinite dimensional manifolcj^ 
with a tangent space T£g = £gx£s and cotangent space T*£j = £j x £^. If F, G 
are multilocal functionals on £g, then their first functional derivatives are smooth, 
so 

F^(p), G (1) M 6 £'s n (D/{u 6 D| (u, ip) = 0 , e £ s }) 

for all <p £ £<j, i.e. ip £ kerP. Since P is a normally hyperbolic operator, one 
can show (see for example m for the proof in a more general setting) that {u £ 

D\ (u,ip) = 0 ,X/p £ £ 5 } = PD, so functional derivatives of multilocal functionals 
are one forms in T(£s x D/PD) C r(T*£s). 

The causal propagator A g induces a Poisson structure on T, which is also well 
defined on the quotient T/To, as To is a Poisson ideal with respect to this structure. 

We can also use A 5 to map one-forms in P(£s x D/PD) C P(T*£s) to one-vectors 
in r(£s x £-s,sc) C r(T£g), where “sc” indicates spacelike-compact support. To 
see how it works, note that A 5 induces an operator D —> £ and ker A 5 = PD, so 
Ag is well defined on equivalence classes in D/PD. To show that As is invertible 
on this space, it remains to show that it is surjective. We recall here the standard 
argument, which can also be found in m■ Let / be a solution with a spacelike- 
compact support, x £ £, and Si, S 2 be Cauchy surfaces such that Si 0 J+(S 2 ) = 0 . 

1 There is another natural way to introduce a smooth manifold structure on £g. We define the atlas 
where charts are given by maps ip + T) —>• £s, sc , p + p Asp, with p £ £s, where £s lSC is the space of 
solutions with compactly supported Cauchy data. We have As : D —> £s jSC and we equip £s. sc with the 
final topology with respect to all curves of the form A 1 —¥ p + As(y>(A)), where A 1 —> p(X) is a smooth 
curve in T>. This gives £s the structure of an affine manifold in the sense of convenient calculus [35] 
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Assume x( x ) = 0 for x £ J_(Ei) and y(:r) = 1 for x £ J+(S 2 ). Then Pxf = 0 
outside of the time slice bounded by Si and S 2 (x =const. there) which implies 
that Pxf has compact support. Hence, 

A S Pxf = A%P X f + A$P (1 - X )f = f- 

We can now assign to a form the vector (A sF^\ •)• On £s jSC we have the 
natural symplectic structure ay: 

(f,g)= [ (f A*dg-*df Ag) = [ (f(d n g)-(d n f)g)dvolv, 

J S Je 

where d n is the normal derivative on S ( d n f = n^d^f, = 0 for £ £ T S, 

n^rifj. = 1.). Obviously, ay extends to a constant 2-form on £g. The relation 
between the 2 -form ay and the “bi-vector field’ll As is given by (see for example m 
for a proof based on the ideas of m 

a 1 (A s F {1 \0 = (F (1 \f) , 

where £ £ T(£s x £s,sc)- In this sense we can think of As as the “inverse” of the 
symplectic structure ay. Setting £ = AgG^ for G £ 3y we obtain 

^(AsFW, A S G (1) ) = (f«,AsG«) ■ 

If S is not quadratic, the situation is more complicated, since S' induces non¬ 
linear equations of motion. It turns out that for many classes of physically interesting 
systems solutions of S'(ip) = 0 develop singularities after a finite time, despite 
starting from smooth Cauchy data. Therefore the space of globally smooth solutions 
£5 might be very small and it does not necessarily capture all the interesting features 
of the theory. Moreover, in general it is not clear if £5 can be equipped with a 
manifold structure in the sense of infinite dimensional differential geometry (see for 
example [2] for the results on the space of solutions of Einstein’s equations). A more 
general structure like a stratified space might be necessary. 

For non-linear equations of motion it is therefore more convenient to replace the 
space of functionals on the space of solutions with the quotient Ts := 3 r /{ (S\ X ), X £ 

V}, where V C T(T£) is the space of vector fields that are derivations of T and the 
duality denoted by (.,.) is the contraction of a 1-form S' £ T(T*£(M)) with a vec¬ 
tor field X. We say that we take the quotient of 3" by the ideal generated by the 
equations of motion. If the equations of motion are linear and normally hyperbolic, 
this ideal coincides with To, so Ts is exactly the space of multilocal functionals on 
£ 5 . In general, our point of view is more in line with the quantum theory and it 
avoids complications related to characterization of the geometrical structure of £s- 
It is also close in spirit to the way one studies varieties in algebraic geometry. 

3.5 Deformation quantization 

Deformation quantization is a method to construct quantum theories from the clas¬ 
sical ones by deforming the commutative product on the space T of functionals to a 

2 Since As is a bi-distribution rather than a smooth function, the map tp 1—As doesn’t induce an 
actual bi-vector field on £, but belongs to a suitable completion of T(A 2 T£). 


13 



non-commutative product * on T[[ft]] (the space of formal power series in h. 


F*G=Y,K nB n( F ,G), ( 40 ) 

n =0 


and we require that 


Bo(F, G) = F ■ G, 

B 1 (F,G)-B 1 (G,F) = ih{F,G}, 

where (F ■ G)(tp) = F(ip)G(tp) is the pointwise product of functionals, {.,.} is the 
Peierls bracket and the second condition is a realization of the idea that in the 
quantum theory one “replaces canonical brackets with commutators”. The existence 
of higher order terms is necessary to avoid the Groenewald-van Hove no-go theorem. 
This result, established first for finite dimensional phase spaces, states that a Dirac 
type quantization prescription is not possible in the strict sense [25, SI]- More 
concretely (see [45]), consider the Lie algebra f) spanned by the canonical coordinate 
and momenta functions q 1 ,... ,q N , Pi, .. . ,pn and 1, equipped with the canonical 
Poisson bracket {.,.} ca n- This is a Lie subalgebra of g = (Pol(T*K JV ), {., .} can ) 
(polynomials on the phase space). The Groenewald-van Hove Theorem states that 
there exists no faithful irreducible representation of \) by operators on a dense domain 
of some Hilbert space which can be extended to a representation of g, so there is no 
quantization map Q from g to the space of operators on some Hilbert space Iff, such 
that 

[Q{f),Q(g)]=ifiQ({f,g}). 

Deformation quantization [10 provides a way out since it weakens the above con¬ 
dition to 

[Q(f), Qia)} = Qdf, 5 ]*) = ifiQiif, g}) + o(h 2 ). 

In field theory, as we have seen in the previous section, the space of functions on 
the IV-dimensional phase space is replaced by T, the space of multilocal functionals 
on £, which is now infinite dimensional. The Poisson structure is provided by the 
Peierls bracket, defined with the use of the causal propagator A 5 . In the simplest 
case, when S is quadratic, one can construct the ^-product using a Moyal-type 
formula. To avoid the functional analytic problems, for the moment we consider 
only regular functionals T r eg, i.e. those for which F t ' n ' > (p) is a smooth compactly 
supported section for all n £ N, <p £ £. On 3yeg[[fi]] we can now define 

OO t]Tl 

(F * G)(ip) = ( Fin) (v), (|A s ) 0n GWM) , (41) 

n =0 


ih / s z \ 

which can be formally written as e 2 \ / F(<p)G(<p') \ v ,'- ip . 

Let us consider the example of a free scalar field and regular functionals of the 
form 

Ff(ip ) = / f(x)tp(x)d 4 x = / f^pdpL, where / £ D. 

J M J 
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We can now define W (/) = exp(iFf) and check that 

((w(/))« (<p), h) = ^ mm*+m u = 0 = fw+^ 

= (i J .fhdfig^ W(/)0p). 


and thus 

((W(/)) W M/") = (ij f h d/j.g'j W (/)(*>). 

Inserting this into the *-product formula, we find, 

W(/)*W(/) = ^ A s (x,y)f(y)f(x)dn g {x)dHg{y) S j W(/+/) 

= e"T A s(/./)w (/ + /), (42) 

which reproduces the Weyl relations. 

Having the interacting theory in mind, we will need to extend the star product 
to functionals more singular than the elements of TVeg, including, in particular, the 
non-linear local functionals. To understand possible obstructions to this extension 
we have to analyze the singularity structure of A 5 . Using the theorem of propagation 
of singularities (see section 0 ), we find that m 

WF(A S ) = {(a -,k-,x,-k') eT*M 2 \(x,k) ~ (x', k ')}, 

where T denotes the tangent bundle minus the zero section and (x,k) ~ ( x',k ') 
means that there exists a lighlike geodesic connecting x and x', to which k is co¬ 
tangent and k 1 is a parallel transport of k. We observe that the WF set of A $ is 
composed of two parts: one with k E {V +) x and another with k E {V-) x , where 
V± is (the dual of) the closed future/past lightcone This observation allows one 
to decompose As into two distributions with WF sets corresponding to these two 
components. Such a decomposition is a local version of the decomposition according 
to positive and negative energies [39]. Therefore we can split As into 

| As = Ait — H , 

where the WF set of Ag is 

WF(A+) = {(*, fc; x, -k’) E T*M 2 \(x, k ) - (x', k 1 ), k £ ( V + ) x } , (43) 

and we also require that As = 2Im(A)t) and that Aj is a distributional bisolution 
to the field equation and is of positive type (i.e. (Ag, / ® /) > 0). On Minkowski 
space one could choose Ait as the Wightman 2-point-function. On general globally 
hyperbolic spacetimes such a decomposition always exists but is not unique. If H 
and H’ correspond to two such choices of decomposition, then H — H’ is a smooth 
symmetric bisolution to the field equations. 

We can now replace |A 5 with Ait in m and the new product, denoted by *h 
can be extended from T reg to defined as the space of functionals with functional 
derivatives satisfying 

WF(FWM)cH n , Vn E N, VV> £ £ , (44) 
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where S„ is an open cone defined as 

S n = T*M n \ {(n,... ,x n ; k u ..., k n )\(k u ...,k n )G (VI uF!) (xi ,..., Xn) } , (45) 

where (V±) x is the closed future/past lightcone understood as a conic subset of 

T* M. 

On IFreg the two star products * and are isomorphic structures and the inter¬ 
twining map is given by 

a H = , (46) 

so that 

Fk H G = a H ((ajf 1 F) * (a^G)) , F,G € T reg . (47) 

In the language of formal deformation quantization one says that products * and 
-kn are related by a gauge transformation, so they provide the same deformation 
quantization. In general a gauge transformation between star products is given by 
f 4 f + where each D n is a differential operator. In our case, 

D n = h.{\( H - H ')iw) n - 

Physically, the transition between * and corresponds to normal ordering, 
so introducing the *#-product is just an algebraic version of Wick’s theorem. As 
stated before, the codomain of an ■ 3>eg —> 3>eg can be “completed” (with the use 
of the Hormander topology El EH) to a larger space 3^ c and we can also build 
a corresponding (sequential) completion a' H 1 (3 r IJC ) of the domain. This amounts 
to extending Tj-eg with all elements of the form linin^oo a] i 1 (F n ) : where ( F n ) is a 
convergent sequence in 3^ c . The quantum algebra 21 of the free theory is defined 
as the space of families Fh, labeled by possible choices of H , where Fh € 21 h = 
(3>c[[fi]])*J?) fulfill the relations 


Fh 1 = olh'-hFr , 


and the product is 

(F k G)h = F h k H Gh- 

We can summarize the relations between the algebraic structures we have introduced 
so far by means of the following diagram: 

(T reg ,*) (Freg,**) 

dense | H dense | PI 

a (T„ c ,k H ) 

A family of coherent states on 2t is obtained by the prescription 

uh,<p(F) = a H (F)(<p) = F H (ip ), 

where ip € £g. This makes sense since Fh is a functional in 3 r /JC , so evaluation at a 
field configuration ip is well defined. 

As an example we can consider the free scalar field with the generalized La- 
grangian 

Lo(f)(<p) = I; J (d^d'V - mV)/ d(j. (48) 

M 
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Let us define 21 as the subalgebra of 21 generated by the Weyl generators W(/) = 
ex~p(iFf). Since 

we conclude that 

SO 

a H (W(/)) = e-t ff ^W(/). 

We can now consider a state obtained be evaluation at tp = 0. We see that 

uh,o CW(/)) = e ~i H UJ) 5 

so H plays the role of the covariance of the state u>h, o- 

Going on-shell corresponds to taking the quotient of 21 by the ideal 2lo generated 
by the elements 

W((n + m 2 )/)-l , /eD . (49) 

Note that S' 0 (tp) = (□ + m 2 )ip , and using partial integration, we can conclude that 
F(jj + m 2 )f( l p) = J S'o (v 3 )/^ = (Sq,/), so taking the quotient by 2to implements the 
free field dynamics. We denote 2t/2to by 21 s 0 j and we see that uih.o is well defined 
on 21 s 0 , as FI is a bisolution for the operator P = □ + in 2 . 

3.6 Interpretation in terms of Kahler geometry 

There is an elegant geometrical interpretation of the structures introduced in the 
previous section. Analogous to Kahler geometry, F[ plays the role of the Riemannian 
metric on y = T>/PT> and the 2-point function = |A s + H is a Hermitian 2-form 

on y. 

The pair (H,As) induces an anti-involution J on y (i.e. J 2 = —1) and if Ait 
is a 2-point function of a quasi-free pure Hadamard state, the triple (H , As, J) is a 
Kahler structure on y. To see how this come about, let us recall some well known 
results (see for example mm for proofs). Let y c denote the complexification of y. 
If A^ and H c are canonical extensions of As and FI to y c , then the following are 
equivalent: 

1. H c + |A§ > 0 on y c , 

2. | (/i, As/ 2 ) | <2(f 1 ,Hf 1 ) 1/2 (f 2 ,Hf 2 ) 1/2 , Sufi G y- 

We can complete y with the product (•>-)h = (■ , iJ.) to a real Hilbert space J~C and 
the inequality 2 implies that As is a bilinear form on K with norm less or equal 2. 
Therefore, there exists an operator A £ 23(J£) with ||A|| < 1 such that 

(f 1 ,A s f 2 ) = 2(f 1 ,Af 2 ) H 

If A has a trivial kernel, then we can just construct the polar decomposition A = 
— J|A| and J satisfies J 2 = —1, so we can use it to equip y with an almost-complex 
structure. More generally, following the proof of theorem 17.12 of m , we can 
define y sg = ker A and y reg = y^. We set A reg = A |"y reg and construct the polar 
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decomposition A reg = — </ reg |A reg |. If the dimension of y sg is even or infinite (which 
is the case in the situation we are interested in), then there exist an orthogonal 
anti-involution J sg on y sg and we set J = J reg © J sg . 

Note that J induces also a complex structure on £g )SC if we set j'Ag/ = Ag,//, 
where / e y. We dehne the holomorphic and anti-holomorphic subspaces of y c as 

Z = {(f-iJf)lfey}, 
z = {(f + iJf)ifey}, 

respectively. Projections onto these subspaces are given by l z = ^(1 — iJ c ) and 
%=i(l+iJ C ) 

The CCR algebra corresponding to (y, Ag) is just the algebra 2lg 0 introduced 
at the end of the previous section. Note that u>h ,o is a state on 21 g with covariance 
H. This state is pure if and only if the triple ( H , Ag, J) is a Kahler structure on y, 
i.e. all three structures are compatible and Ag o J = 2 H. We can now decompose 
As in the holomorphic basis. A straightforward computation shows that 

<%/ 1 ,A+(l z / 2 )) = (/ 1 ,A+/ 2 ) , 

where /i, / 2 £ y c and remaining components vanish, so in the holomorphic basis 
Ag is represented by 

0 0 \ 

Ag Oj ’ 

so it acts only on the holomorphic part of the first argument and the anti-holomorphic 
part of the second argument. 


4 Time ordered products, and the perturbative con¬ 
struction of local nets 

In the previous section we were concerned only with the quantization of free theories 
(quadratic actions). Given an arbitrary action S we first split S = Sq + Si , where So 
is quadratic. We already know how to quantize the classical model defined by So, so 
now is the time to introduce the interaction. We will do it in this section, following 
the ideas of H E123 M M HQ, but before we start, we give a heuristic argument 
justifying our construction. The idea is to use the analogy with the interaction 
picture of quantum mechanics. Let Hq be the Hamiltonian operator of the free 
theory and let H t j = — f K :£/(0, x): da t be the interaction Hamiltonian, where :£/: 
is the normal-ordered Lagrangian density, constructed from the classical quantity 
Lj and K is some compact subset of E (as explained in section ED- The rigorous 
“smoothed-out” version of the Hamiltonian quantization will be given in section [6] 
We would like to use the Dyson formula and define the time evolution operator 
as a time ordered exponential, i.e. 

U{t,s) = e ^o e -i(t-s)(ffo+ff/) e -i S ffo = 

OO n 

1 + V'^-r / T(:£/(cci):... :£/(x n ):)d 4rl x, 

„=i n] Amixr 3 )" 
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where 

x i-> £ f (x) = e iH ° x ° :L I (0,x):e- iHoX ° 

is an operator-valued function and T denotes time-ordering. Heuristically, one could 
use the unitary map defined above to obtain interacting fields as 

ipi(x) = U(x°, s)^ 1 <p(x)U(x °, s) = U(t, s)~ 1 U{t, x°)ip(x)U(x°, s) , (50) 

where s < x° < t. 

There are, however, serious problems with this heuristic formula. Firstly, typical 
Lagrangian densities, e.g. :£/(x): = :</?(x) 4 :, can not be restricted to So as operator 
valued distributions. This is the source of the so called UV problem. Moreover, as 
mentioned before, having the sharp cutoff function in the Lagrangian and Hamilto¬ 
nian (like in (1 25 [) 1 leads to additional divergences (Stiickelberg divergences). Finally 
there is the adiabatic limit problem related to the fact that the integral over x does 
not exist. Last but not least, the overall sum might not converge. 

4.1 Causal perturbation theory 

Some of these problems mentioned in the introduction can be easily dealt with by 
a slight modification of the above ansatz. For example, we avoid the Stiickelberg 
divergences by replacing the sharp cutoffs with smooth test functions. The UV 
problem is solved by using causal perturbation theory in the sense of Epstein and 
Glaser [HsS| . In this method one switches the interaction on only in a compact region 
of spacetime and then takes the adiabatic limit (understood as a certain inductive 
limit) on the level of interacting observable algebras. These modifications of the 
Dyson formula ansatz lead to the definition of the formal S-matrix: 

oo , n 

S<»> = 1 + ES 

n—1 

where g is a test density. In order to make this formula well defined, we need to make 
sense of the time-ordered products of :Lj(xi):. This will be done by Epstein-Glaser 
renormalization. Finally, the formula (1501) has to be reinterpreted as a definition of 
a distribution, rather than a function. Hence, for a test density / we obtain 

OO .yi ft 

f( x ) ipi (x) = S(g)- 1 J2 l - / f(x)g(x 1 )...g(x n )Ty(x)£ I (x 1 )...£i(x n ) = 

n=0 n ' 

= -^ s (g)~ ls (g,Xf)\ x=0 , 

where S(g, f) is the formal S-matrix with the Lagrangian density gLi + ftp. This is 
the so called Bogoliubov’s formula [7]. 

We are now left with the problem of defining the time-ordered products on 
Wick-ordered quantities :£/(x):. We have already mentioned in section 1531 that 
the normal ordering corresponds to passing between the star product on 21 h 
and * on 21. Note that elements of 21 h are functionals on £s, so we can identify 
classical quantities in Ti oc with quantum ones by means of 7^ : Ti oc —>• 21 h defined 
by 7i = id. Composing with afj 1 we obtain a map Ti : $ —> 21, 7± = afj 1 o 


J g(x i)... g{x n )T(\£ I {x i):... :L I (x n ):) 
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which maps “classical” to “quantum”. This map is interpreted as the normal ordering 
and we can now make an identification 


:F: = Ti F, F e T loc . 


In the context of local covariance, this choice of normal ordering is not the most 
optimal one. This is because a family of Hadamard states cannot be chosen in a 
covariant way (i.e. compatible with embedding of globally hyperbolic spectimes), 
but a family of Hadamard parametrices can. The latter are bi-solutions of the 
linearized equations of motion only up to smooth terms. It is, therefore, more 
appropriate to define the normal ordering by a prescription where only the singular 
part of H is subtracted from the correlation function of two fields, as opposed to 
the prescription where one subtracts the full H. Concretely, we set Tf* = a w 
so Ti = , where w is the smooth part of the Hadamard 2-point function 

(see [M] and m for a recent review). More precisely, this has to be understood 
as linyy^oo F for F € Tioc and this limit makes sense, because the series 

converges after finitely many steps. The function u>n appearing in this prescription 
is 2 N + 1 times continuously differentiable and it appears in the 2 -point function 
as Ag = W^ n9 + wn- The singular part W^ n9 is of the form “^ + v In a'\ with 
<j{x, y ) denoting the square of the length of the geodesic connecting x and y and with 
geometrical determined smooth functions u and v. For a more precise definition of 
what is the Hadamard form for of a 2-point function, see for example [34] or a recent 
review m • 

More concretely, for a density of the form 

/ d n 

f{x)—A(x)(v>)\ x= 0 d 4 x, 


where A(x)(tp) = (here p is a polynomial in covariant derivatives) we 

define 

/ J7l 

f^ A n\ x= 0 d 4 x, 

where 

A h (x) = e ^ v) ® piv)wN{x ' x) A{x ), 

Unfortunately, the modifications which we have done so far do not render the 
time-ordered products well defined. Heuristically, we would like the time-ordered 
product of two functionals to be 0, i.e. 

(F • t G) = eK^' A i > o^)F(^)G(^ , )l^= v , 

where A<? = ^(A<j o + A^ o ) is the Dirac propagator. This makes sense if both F and 
G are regular functionals (i.e. elements of 3>eg)- This indeed provides the correct 
notion of time-ordering, since 


( F *G if supp G -< supp F, 
( G * F if supp F supp G , 


(51) 


where the relation means “not later than” i.e. there exists a Cauchy surface 
which separates supp G and supp F and in the first case supp F is in the future of 
this surface and in the second case it’s in the past. 
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The time ordered product defined by ED is associative, commutative and iso¬ 
morphic to the point-wise product by means of 



f.iG = T(r 1 F.r 1 G) , 

(52) 

where 




ih(A% ,AL\ 

T = e s o’JJz' 

(53) 

or more precisely 




OO i— yi 

( xf)m = y, ^ {(iAgr n ,F^\v)) ■ 

n— 0 

The linear operator T defined above is sometimes called the “time-ordering operator” 
and it is interpreted as a map which goes from the “classical” to the “quantum”, i.e. 

(Treg,-) ^ (a l-eg) "?) 

classical quantum ’ 

where 2t reg C a is the range of T. Note that on the quantum side we have two 
products. Using the time-ordered product we can express the formal S'-matrix S : 
3 reg[[fi]] —i► 3 reg[[fi]] as the time ordered exponential: 

S(V) = ef /R = T(e T_liy / fi ) . (54) 

According to our interpretation of T, § is a map on the “quantum” algebra a to 
itself. Interacting fields are obtained by means of the Bogoliubov formula, which 
reads 

R V (F) = (S(V )*- 1 *S(U + AF))| a=0 = 

= (e’7 /R )* _1 * (e l 7 /n F) . (55) 

We interpret Ry{F) as the interacting quantity corresponding to F. We can also 
define the interacting star product as 

F -k V G = Ry\R v F * R V G). 

The interacting theory is given in terms of the algebra (T reg ,*y) and Ry acts as 
the intertwining map between the free quantum theory and the interacting quantum 
theory, i.e. 

(~Ueg; ■) t (a reg ,*,’ t) r v 1 (a reg ,*y) , . 

classical fre ) e interacting ' l > 

quantum quantum 

All these formulas make sense if we restrict ourselves to regular functionals. This 
is, however, not satisfactory for our purposes, since typical interactions are local and 
non-linear, hence not regular. In the first attempt we could try to pass to a different 
star product, which amounts to replacing A s 0 by Ag g and At? by the Feynman 
propagator Ag g = iAg g + Ft , so our diagram gets modified to 

((Jreg, - ) T H ' \ a H 1 , (^4eg)*)' T ) Ry 1 ^ ( 2 lreg>*v) 

classical 1 reg ’ Hr,jH) interacting* 

quantum quantum 
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where 7 H = e lh< ' As °’w T \ so 7 = a T 1 o 7 H . This modification of the formalism, 
however, doesn’t solve the problem yet. To extend our formalism to arbitrary local 
functionals, we need to perform the renormalization. The difficulty which we have 
to face is the fact that the WF set of Ag 0 at 0 is like the WF set of the Dirac delta 
and therefore the tensor powers of Ag cannot be contracted with derivatives of 
local functionals. 

However, there is a way to extend 7 H to local functionals. First we extent 7 H to 
i?ioc by setting 7 H = Tf. We discuss here only the Minkowski spacetime situation, 
so we can set Tf = id. The subspace 7 H (7\ oc ) C 21 will be denoted by 2lf c . Let us 
define the n-th order time-ordered product as 

(Fi, • • •, F n ) = F V7h ... -t H F n , 

whenever it exists. It is well defined for F\,... ,F n £ Ti oc with pairwise disjoint 
supports and we will denote this domain of definition by (Ti oc )ff s . Moreover 

T- (F u ...,F n ) = 7Z(F 1 ,..., F k ) * H 7 f_ fc (F fc+1 , ...,F n ), (T 1) 

if the supports suppFi, i = 1 ,... ,k of the first k entries do not intersect the past 
of the supports supp Fj, j = k + 1,..., n of the last n — k entries. This property is 
called the causal factorisation property. We will take it as an axiom that we want 
to impose while extending time-ordered products to arbitrary local arguments. The 
other axioms include 

(T 2) Starting element: T- = 1, Tf = id, 

(T 3) Symmetry: Each Tf is symmetric (graded symmetric if Fermions are present). 

(T 4) p -Locality: Tf (F ±,..., F n ), as a functional on £, depends on ip only via the 
functional derivatives of F \,..., F n . 

In the seminal paper [18] , Epstein and Glaser have shown that such a family of maps 
exists and non-uniqueness in defining Tf’s is fully absorbed into adding multilinear 
maps Z n : 2lj®f ->■ 2li oc , he. 

T^fFi, ...,F n ) = Tf (F u ...,F n ) + Z n (Fi, ...,F n ), 

where {Tf } n gN an d {T H n }„ g N are two choice of time-ordered products that coincide 
up to order n — 1. The renormalized S-matrix is now defined by 

oo oo 

s (^) = E ■ ■ • < *0 = E ° Vn (* H V,..., auV). 

n =0 n —0 

The causal factorisation property for time ordered products implies that the S-matrix 
satisfies Bogoliubov’s factorization relation 

S(Vi + V 2 + V 3 ) = §(Vi + V 2 )S{V 2 )- 1 S(V 2 + V 3 ) (57) 

if the support of V\ does not intersect the past of the support of V 3 . 

We can also define the renormalized map T : T —» 21 by T = 0 ?l aff o Tf o m , 
where m~ l : 7 S'7 is the inverse of the multiplication, as defined in [20] 

and T® is the space of local functionals that vanish at 0. The renormalized time 
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ordered product -t is now a binary operation defined on the domain D<j = ‘T(T). 

Analogously to the diagram (1561) . we obtain now 

(3V) T (%)*! "*) Ry 1 (21,*V) /-O') 

classical fr< ; e interacting > 1 ) 

quantum quantum 

with the caveat that -t is well defined on D-j C 21. 

We will now discuss in detail the ambiguity arising in defining T n ’s. In physics 
this is known as the renormalization ambiguity. To understand it better and to 
relate it with the notion of the renormalization group, we first define a map 2 : 

2lioc[[fi]] —t 2ti oc [[ft]] by summing up all the Z n ’s relating two chosen prescriptions 
to define the time-ordered products. For any two choices of T n ’s the corresponding 
map Z has the following properties: 

(Z 1) 2,(0) = 0, 

(Z 2) 2«(0) = id, 

(Z 3) 2 = id + 0(H), 

(Z 4) 2 (F + G + H) = 2 (F + G) + Z(G + H) - 2(G), if supp F n supp G, 

(Z 5) f = 0. 

The group of formal diffeomorphisms of 2ti oc [[/i]] that fulfill (Z [TJ) - (Z [5]) is called 
the Stiickelberg-Petermann renormalization group Ik. There is a relation between 
the formal S-matrices and elements of 1R provided by the main theorem of renormal¬ 
ization mm. It states that for two S-matrices § and S, built from time ordered 
products satisfying the axioms (T 1) - (T [4}, there exists 2 £ !k such that 

S = So2, (59) 

where 2 £ 1R and conversely, if § is an S-matrix satisfying the axioms (T 1) - (T[4l) 
and Z £ Jl then also § fulfills the axioms. 

4.2 Methods for explicit construction of time-ordered prod¬ 
ucts 

The proof of existence of time-ordered products with properties (T 1) - (T 0) 
given in )l 8 j is rather abstract and relies on an inductive argument. For practical 
purposes an existence result is not sufficient and one would like to obtain some 
explicit formulas for T ra ’s. In this section we will review results which show that 
the problem of constructing time-ordered products reduces to extending certain 
distributions. Subsequently, we will give some concrete computational prescriptions 
for constructing such extensions. 

We start with an example. Let F = \ f p 2 fdp, G = ~ f ip 2 gdp<, f,g £ D. If 
supp <7 fl supp / = 0 , then the time ordered product -<j of F and G is given by 

T 2 (F,G)M = 

= (F"jG)(ip) = F{ip)G(tp)+ih J <p(x)ip(y)f(x)g{y)A% a (x,y)dp,(x)d[i(y)-^y j A% 0 (x, y) 2 f(x)g{y)dy(x)dp(y) 

In the least term of the expression above we have a pointwise product of a distribu¬ 
tion with itself. This could potentially cause problems. If x ^ y, then if (x,k) and 
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(y, —k') belong to the wave front set, then k, —k! are cotangent to a null geodesics 
connecting x and y. Moreover, k is future directed if x is in the future of y and 
past directed otherwise, so the sum of two such covectors doesn’t vanish. Hence, 
the condition on the multiplicability of distributions presented in section [ 2 ] implies 
that (Ag o ) 2 as a distribution is well defined on the complement of the diagonal 
{(x,x)\x £ M}. Let us now consider what happens on the diagonal. There, the 
only restriction is k — —k', hence the sum of WF(Ag o ) with itself contains the zero 
section of the cotangent bundle at the diagonal. The problem of defining Tf (F, G) 
reduces now to the problem of extension of Ag o to a distribution defined everywhere. 

This generalizes, and the construction of T"’s reduces to extending numerical 
distributions defined everywhere outside certain subdiagonals in M n . The construc¬ 
tion proceeds recursively and, having constructed the time-ordered products of order 
k < n, at order n one is left with the problem of extending a distribution defined 
everywhere outside the thin diagonal of M n . On Minkowski spacetime, exploiting 
the translational symmetry of M, this reduces to extending a numerical distribution 
defined everywhere outside 0. One way of constructing explicitly such distributional 
extensions relies on the so called splitting method (see for example HO])- Here we 
will take a different approach, based on the notion of Steinmann’s scaling degree 
HU. Here is the definition: 

Definition 1. Let U C R n be a scale invariant open subset (i.e. XU = U for 
A > 0), and let t £ D'(U) be a distribution on U. Let t\(x) = t(Xx) be the scaled 
distribution. The scaling degree sd of t is 

sdt = inf {<5 £ Ml lim X s t\ = 0} . (60) 

A—>o 

The degree of divergence, another important concept used often in the literature, 
is defined as: 

div(f) = sd(i) — n. 

The crucial result which allows us to construct time-ordered products is stated in 
the following theorem: 

Theorem 1. Let t £ D(R” \ {0}) with scaling degree sdt < 00 . Then there exists 
an extension oft to an everywhere defined distribution with the same scaling degree. 
The extension is unique up to the addition of a derivative P(d)8 of the delta function, 
where P is a polynomial with degree bounded by div(t) (hence vanishes for sd t < n). 

In the example presented at the beginning of this subsection, the scaling degree of 
in 4 dimensions is 4, so the extension exists and is unique up to the addition 
of a multiple of the delta function. 

The result above allows in principle to extend all the numerical distributions we 
need for the construction of time-ordered products. However, the computations can 
in general get very complicated, so it is convenient to formulate the combinatorics 
underlying our construction in terms of Feynman graphs. In the pAQFT framework, 
these are not fundamental objects, but instead they are derived (together with the 
corresponding Feynman rules) from time-ordered products. 

Time-ordered products should be maps from T®" to 3>ic[[ft]] and, as indicated 
in the previous section, they are obtained by extending non-renormalized expressions 
that are originally defined only on (Tioc)^- Let us consider F = F\ ® • • • (8) F n £ 
(^loc)®^ with the corresponding Wick-ordered quantities are elements of 2li oc given 
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by A\ = 7Fi,... ,A n = T F n £ ^loc- Note that F induces a map from £ ra to R 
by F(ipi ,..., ) = -Fi(^i)' •' F n (ip n ). When we talk about functionals on £ we will 

denote the variable by ip and for functionals on £™ we take an n-tuple {p\ 7 ip n ). 

Let us denote D %] = ih(Ag o , s 5 s ) and D = ih{ Ag o , ^-y). The Leibniz rule for 
differentiation can be formulated as 


<5 

— o m n = m n o 
op 



(61) 


where m n is the pointwise multiplication of n arguments, or in other words, a pull¬ 
back through the diagonal map £ —> £ n , ip i-> (ip, The Leibniz rule implies 

that the non-renormalized expression for 7 H satisfies 

7 H o m n = e~ o m n = m n o e^’<j ; 


Hence 


-hD i 


Fi-<jh■ ■ ■•<jhF ti = e 2 om„(e 2 




F n ) = m n oe L -<^ ij (Fi,... F n ) = m n oT n (Fi,. 


We can now use an identity 


OO 

-liv 

lij —0 


<j 


n 


i 


(62) 


to express time ordered products in terms of graphs. Let be the set of all graphs 
with vertex set V(T) = {1,.. .n} and Uj the number of lines e £ E(T) connecting 
the vertices i and j. We set = Iji for i > j and la =0. If e connects i and j we 
set de := {i 7 j}- Then 

T n =J2 T r, (63) 

res„ 


where 


with 


and 


T r = 


1 


(5r = 


Sym(T) 

^2 |s(r)| 


(tr, <^r) 


riiey(r) Yie-.iede ^ < fi{ X e,i) 


tr = ]^[ hAp(xe t i,i £ de) 
e£B( r) 


(64) 


(65) 


The, so called, symmetry factor Sym is the number of possible permutations of 
lines joining the same two vertices, Sym(T) = . lij\. Note that Tr is a map from 

(l^ioc)®^ to C°°(£l v ^,R)[[fi]], where <8 >V means that the factors in the tensor product 
are numbered by vertices and to a vertex v £ V (T) we assign the variable ip v . The 
renormalization problem is now the problem to extend T n ’s to maps on (Tio C )®" 
and this can be achieved by extending all the maps Tr and using formula (1631) . 

First we note that functional derivatives of local functionals are of the form 


F {l) (p)(x 1 ,...,xi) 


r N l 

= m > ■ ■ ■ > d *i ) n 5 ( y 

^ j=1 i= 1 


x i)dp(y ), 


( 66 ) 


■■F n ). 
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where N £ N, p 3 's are polynomials in partial derivatives and g 3 [p\ are (^-dependent 
test functions. The representation above is not unique, since some of the partial 
derivatives d Xi can be replaced with d y and applied to g 3 [tp\. Another representation 
of (ip) is obtained by performing the integral above and using the centre of mass 
and relative coordinates: 

F (l) ... ,xi) = '^2fg[p\{z)d f3 S{x iel ) (67) 

P 


where j3 £ test functions fp[<p\{x) £ D are now (^-dependent functions of 

the center of mass coordinate 2 = {x\ + • • • + Xk)/k and x rel = {x\ — z, ... ,Xk — z) 
denotes the relative coordinates. 

Using (1551) we see that the functional differential operator <5p applied toFe fFj®" 
yields, at any n-tuple of field configurations ..., <p n ), a compactly supported dis¬ 
tribution in the variables x e ,i, i £ de , e £ E(T) with support on the partial diagonal 
Ar = {x e p = Xf,i,i £ de fl df, e, / £ E(T)} C and with a wavefront set 

perpendicular to TAp. Note that the partial diagonal Ap can be parametrized using 
the center of mass coordinates 

v valence(n) e ’ v ’ 

ejvGde 

assigned to each vertex. The remaining relative coordinates are x™l = x e:V — z v , 
where v £ V(T), e £ £?(r) and v £ de. Obviously, we have Y^ e \vede x *e}v = 0 f° r all 
v £ V'(T). In this parametrization SrF can be written as a finite sum 

SrF = J2 /^Aei, 

finite 


where /3 £ Nq V each fP(tpi ,... ,<p n ) is a test function on Ap and <5 re i is the Dirac 
delta distribution in relative coordinates, i.e. S Tf ,i(q) = o(0, ...,0), where q is a 
function of £ V(T),e £ E(T)). 

We can simplify our notation even further. Let Yp denote the vector space 
spanned by derivatives of the Dirac delta distributions dpS le 1 , where /3 £ Nq 1 
Obviously, Yp is graded by |/3|. Let CD(Ap,Yp) denote the graded space of test 
functions on Ap with values in Yp. With this notation we have SrF £ D(Ap,Yp) 
and if F £ (Tioc)®^, then SrF is supported on Ap \DIAG, where DIAG is the large 
diagonal: 

DIAG = {z £ Ap| 3v, w £ V (T), v ^ w : z v = z w } . 

We can now write (1551) in the form 


1 

Sym(T) 


{tr,Sr) = ^2 (f^dptirehtr) 

finite 


where tp is now written in terms of centre of mass and relative coordinates. To see 
that this expression is well defined, note that we can move all the partial derivatives 
dp to tp by formal partial integration. Then the contraction with (5 re i is just the 
pullback through the diagonal map map pr : Ap —> M 2 l E l r ll by 


(pr(z)) e ,v = z v if v e de. 
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From the wavefront set properties of Ag o , we deduce that the pullback pj of each 
tp = dptr is a well defined distribution on Ar\DIAG, so ((Ml) makes sense if F £ 
(Jioc)® d " s , as expected. We conclude that £r £ ©'(Ar\DIAG, Yr), where the duality 
between tr and a test function / = ^finite is given by 

0 


The renormalization problem now reduces to finding the extensions of £p, so that 
fp gets extended to an element of ©'(Ar,Fr)- The solution to this problem is 
obtained by using the inductive procedure of Epstein and Glaser. The induction 
step works as follows: if tr' is known for all graphs T' with fewer vertices than 
T, then tr can be uniquely defined for all disconnected , all connected one particle 
reducible and all one particle irreducible one vertex reducible graphs. Graphs which 
are irreducible and do not contain any non-trivial irreducible subgraphs are called 
EG-primitive. For the remaining graphs, called EG-irreducible , tr is defined uniquely 
on all / £ ©(Ar,Fr) of the form above where fp vanishes together with all its 
derivatives of order < up + |/?| on the thin diagonal of Ap. Here 

uj t = (d- 2)|£(r)| -d(|V(r)| - 1) 

is the degree of divergence of the graph T. We denote this subspace by © wp (Ar, hr)- 
Graphs which are irreducible and do not contain any non-trivial irreducible sub¬ 
graphs are called EG-primitive. Renormalization amounts to project a generic / to 
this subspace by a translation invariant projection Wp : ©(Ar, Fr) —t © wr (Ar, Fr). 
Different renormalization schemes differ by different choices of the projections Wr 
(see m for details). 

On Minkowski spacetime we have further simplifications. By exploiting the 
translation invariance we find that, at each step of the recursive construction of 
time-ordered products, the renormalization problem reduces to the problem of ex¬ 
tension of some distribution defined everywhere outside the origin, so this is what 
we will focus on now. 

For concrete computations it is convenient to construct these extensions with 
the use of regularization. Let us first define the notion of a regularization of a 
distribution. Let t £ ©'(R d \ {0}), d £ N, be a distribution with degree of divergence 
w, and by t £ ©(,(R d ) we denote the unique extension of t with the same degree 
of divergence. A family of distributions {^}(ef!\{o}i ^ £ ©'(R"), with Q C C a 
neighborhood of the origin, is called a regularization of t, if 

Mg £ © A (R d ) : lim (t<,g) = (t,g). (68) 

We say that the regularization {f^} is called analytic, if for all functions / £ ©(R n ) 
the map 

H\{0}3C (69) 

is analytic with a pole of finite order at the origin. The regularization {f*’} is called 
finite, if the limit lim^o^) /) £ C exists V/ £ ©(R d ). 

For a finite regularization the limit lim^o A is, as expected, a solution t of the 
extension (renormalization) problem. Given a regularization {t.^} of f, it follows 
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from (RISll that for any projection W : T> —> D^ 

(i,Wf) = lim{t c ,Wf) V/eB(R' n ). (70) 

It was shown in m that any extension t £ 2)'(R d ) of t with the same scaling degree 
is of the form (t, f) = (t, Wf) with some W-projection of the form 

Wf:=f- J2 / (Q) (°)«’«. (71) 

l«l<A 

where w a £ D(R d ) such that for all multiindices /? £ Ng with \j3\ < oj we have 
d^w a ( 0) = 6%, |a| < w Hence 


(t, Wf) = lim 


(t ( J) 


E (i C .^a)/ ( “ } ( 0) 

|o;|<sd(t)— n 


(72) 


In general, we cannot split the limit on the right hand side into two well defined 
terms. However, if the regularization {t^,( £ II \ {0}} is analytic, then we can 
expand each term into a Laurent series around £ = 0, and because the overall limit 
is finite, the principal parts (pp) of these two Laurent series must be the same. This 
means that the principal part of any analytic regularization {t^} of a distribution 
t £ D'(R d \ {0}) is a local distribution of order sd(f) — d. Following [17|, we can now 
give a definition of the minimal subtraction in the EG framework. 

Definition 2 (Minimal Subtraction). The regular part (rp = 1 — pp,) of any analytic 
regularization {t^} of a distribution t £ ©'(R^ \ {0}) defines by 

(t MS ,/):=limrp((ff,/)) (73) 

C->o 


an extension of t with the same scaling degree, sd(f MS ) = sd(t). The extension t MS 
defined by m is called the “minimal subtraction”. 


4.3 Interacting theories 

Let us now discuss the problem of constructing interacting nets of observables. We 
start from a space T) n of functions / : M —>■ R” with compact support. We assume 
that we have unitaries S(f), f £ T) n with S'(O) = 0, which generate a *-subalgebra 
21 of 21 and satisfy for f,g,h£T> Bogoliubov’s factorization relation 

S(f + g + h) = S(f + g)S(g)- 1 S(g + h) 

if the past J_ of supp h does not intersect supp / (or, equivalently, if the future J + of 
supp / does not intersect supp h). We can obtain these as formal S-matrices S(f) = 
§(V(f)), discussed in the previous section (see property (1571) ), for a generalized 
Lagrangian V(f) = a H > where / £ D n and each Aj(x) 

is a local function ip £ £. Typically 4' are polynomial and they represent Lagrangian 
densities of various interaction terms that one can add to the free action Sq. 

We also assume that the translation group of Minkowski space acts by automor¬ 
phisms a x on 21 such that 

ax(S(f)) = S(f x ) , f x {y) = f(y - x) . 
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Obviously, this is also satisfied for the S-matrices discussed so far. Under these 
general assumptions, we define local algebras 21(0), 0 C M, as the *-subalgebras of 
21 generated by S(/),supp/ C 0 and obtain a translation covariant Haag-Kastler 
net on Minkowski space. To justify this claim, we will now check that all the axioms 
are satisfied. 

Isotony and Covariance are obvious, and Locality follows from the fact that for 
functions /, g with spacelike separated supports 

supp / (~l ,J± (supp g) = 0 (74) 

and hence 

S(f)S(g) = S(f + g) = S(g)S(f) . (75) 

The crucial observation is now that the map / H > S(f) induces a large family 
of objects that satisfy Bogoliubov’s factorisation relation, which are labeled by test 
functions g £ CD™, namely the relative S-matrices 

f^S g (f) = S(g)- 1 S(g + f) . 

We can choose Ao(x) = £i(x) to be the Lagrangian density of the interaction term. 
Then, for g = ( 50 , 0 ,... ,0), we obtain V(g) = f £/g 0 d/n = T/(ffo), where g 0 £ CD. 
Note that S(g + A/) = §(a H (Lj(go) + A JT f Ajfjd/j,)), so the derivative of S with 
respect to A is just the retarded field Rx, I ( go )(V(f))- Let us now prove that the 
causal factorisation property indeed holds for S g (f). Let /, h £ CD™ such that supp / 
does not intersect J_(supph). Let g,g' £ CD™. Then 

Sg(f + 9 + h) = S{g) 1 S{f + (g + g) + h) 

= S(g)~ 1 S(f + {g + g'))S(g + g'^Siig + g') + h) 

= S g (f + g')Sg(g')- 1 S(g)~ 1 S(g) S g (g f + h) . □ 

=1 

We consider S g (f) as the retarded observable S(f) under the influence of the 
interaction Li(go). The Haag-Kastler net 2l s of the interacting theory is then de¬ 
fined by the local algebras 2l ff (0) which are generated by the relative S-matrices 
S g (f), supp / C 0. These can indeed be interpreted as retarded observables, as 
S g (f ) depends only on the behavior of g in the past of supp/. More precisely, 
supp(g — g') PI J_(supp /) = 0 implies 

S g {f) = S{g)~ 1 S((g - g') + g' + /) = 

= S{g)~ x S{{g - g') + gfW)- 1 Stf + f) = S g ,(f) . 

The second observation is that S g (f) depends on the behavior of g outside of the 
future of supp / via a (formal) unitary transformation which does not depend on /. 
Namely, supp(g — g') (~l J+(supp /) = 0 implies 

S g (f) = S(gr 1 S(f + g' + (g-g')) = 

= S(g)~ 1 S(f + g')S(g')~ 1 S(g' + (g - g')) = 

= S(g)~ 1 S(g')S(g')~ 1 S(f + g')S g> (g - g') = 

= Ad Sg'(g-g')-\SAf)) . 
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Hence the structure of local algebras depends only locally on the interaction. This 
allows to perform the adiabatic limit directly on the level of local algebras. 

In the next step we want to remove the restriction to interactions with compact 
support. Let G : M —> R” be smooth and 0 be bounded. Set 

[G]o = {g £ T) n \g = G on a neighborhood of J+(0) fl J_(0)} . 

We consider the 21-valued maps 

SG,o(f)--[G\oBg*S g (f)e%. 

The local algebra 21g(0) is defined to be the algebra generated by Sa,o(f ), supp / C 
0. Note that the evaluation maps 

7 gG : S G ,o(f) ^ S g (f) 

extend to isomorphisms of 21g(0) and 2l s (0) for every g £ [G]q. 

The local net is now defined by the embeddings *o 2 Oi f° r Oi C O 2 

«o 2 Oi : SoMf) 5g,o 2 (/) 

for / £ T> n with supp/ C Oi. Let 21 g be the inductive limit with embeddings 

*0 : 21g(0) —> 21 g 


and we set 

S G (f) = i 0 (S G , 0 (f)). 

We are now ready to prove a crucial theorem about the net 0 21 g(0). 

Theorem 2. Let G be translation invariant. Then the net becomes translation 
covariant by setting 

a x (S G (f)) = S G (f x ) . 

Proof. We have to prove that a x extends to an isomorphism from 21 g(0) —5- 21g(0 + 
x). Let Oi D 0 U 0 — x and g £ [G]oi- Then g,g x £ [G]o and g x = g + h+ + h- 
with supp/iff. n J T (0) = 0. By causal factorization 

a x = % G 0 AdU g (x) oa x o 7 9 g 

with U g (x) =S g (h x _). □ □ 

5 Time-slice axiom, operator product expansions, 
and the renormalization group 

We have seen that, starting from a free QFT and a definition of a time ordered 
product satisfying the axioms of Section 14.31 we can construct a local net (in the 
sense of formal power series) satisfying the Haag-Kastler axioms of Isotony, Locality 
and Covariance. In this section we want to analyze the net in more detail. 

First we investigate whether the net satisfies the time-slice axiom. This can be 
done for the case that the net is defined on a generic Lorentzian globally hyperbolic 
spacetime M. It is known since a long time [22] that the free theory generated by 
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linear functionals, modulo the ideal of the free field equation, satisfies this axiom, 
and by using the techniques of microlocal analysis, this result can be extended to 
the net 21 generated by elements of the form ct B F, where F £ 5^ is a microcausal 
functional E3 HU- In eu it was shown that this implies that also the net 21 g 
introduced in the previous section satisfies the axiom. The argument relies on the 
fact that the algebra of the interacting theory associated to some bounded region 
can be constructed as a subalgebra of the free theory for a slightly larger region, 
and vice versa. 

The problem is that these subalgebras are fixed only up to unitary equivalence, so 
one has in addition to show that these unitary transformations can be appropriately 
fixed. We use the fact that the relative S-matrices S g (f) are well defined also for 
test functions g with non-compact support provided the support is past compact, 
i.e. supp g fl J_( x) is compact for all x £ M. 

Let E be a Cauchy surface of M and N a neighborhood of E. Let OcMbe 
relatively compact. We choose a Cauchy surface E_ such that OUTV C J+(E_) and 
a smooth function \ with past compact support such that supp(l — x) C J_(E_). 
We want to prove that 21 gx(Q) C 21 Gx (N). 

By construction of the interacting theory we see immediately that 21 Gx (0) C 
21 (M) holds. Due to the time slice property of the free theory, 21 (M) = 21 (IV 7 ) 
for each neighborhood N' of E. We now construct within the algebra 21 gx(A 0 an 
algebra which is isomorphic to 21 (IV 7 ) for a sufficiently small E C N' C N. For this 
purpose we choose another smooth function x' with support contained in J + (iV) 
and with supp(l — x') C J-(N'). Let now supp / C N'. Then the unitaries 

S G ( x -x')(f) = S G X W)~ ls GxW + /) > with g' = G\' on J_(supp/),supp g' C N , 

(76) 

generate an algebra isomorphic to 21 (IV') within %l Gx (N). The map 

a : S(f) -»■ S G{x _ xl) (f) = Ad(S(g - g'))-\S(f)) (77) 

with g = G\ on J_(supp/ U suppg') extends to an injective homomorphism from 
21 (M) into 2l Gx (N). Since 21 Gx {N) C 21 (M), a is an endomorphism of 21 (M). We 
show that it is even an automorphism. For this purpose we construct the inverse of 
a. By exploiting the time slice property of the free theory, we can restrict ourselves 
to elements S(f) with supp/ C J_(E 0 ). On these elements we have 

a" 1 (5(/)) = Ad (S(g - g'))(S(f)) = S(g - g'+ f)S(g - g')~ 1 (78) 

where g — g' = G(x ~ x') on T+(supp /). We conclude that 21 Gx {N) = 21 (M). This 
proves the claim. 

Another general property of the interacting net is the existence of an operator 
product expansion |28| . In the case of the product of two fields A and B it is an 
expansion 

A(x)B(y) ~ ^2 Cab(x, y)yk{x) (79) 

k 

with distributions C\ B and a basis of local fields ipk, ordered with respect to the 
scaling dimension. This is an asymptotic expansion in the sense that after evaluation 
in a state coming from a Hadamard state of the free field, the difference between the 
right hand side of the relation and the left hand side, truncated at some k 7 tends to 
zero as x —> y, with an order depending on k. 
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The third property we look at is the behavior of the theory at different scales. 
In the standard formalism of QFT, one formulates this as a property of vacuum 
expectation values of products or time ordered products of fields, or one uses the 
concept of the so-called effective action. In this formulation one has to have control 
over the existence and uniqueness of the vacuum state. In the algebraic approach 
one can instead derive a relation between local nets. Namely given a local net 
0 i—^ 2li(0) one obtains another net by scaling the regions, 

21a( 0) = Sir (AO) . (80) 

If the net depends on some parameters ( m,g ), one can compensate the scaling by 
changing the parameters. One obtains the algebraic Callan-Symanzik equation [9] 

21 ™’ 9 ^ 21^(0,30) (81) 

The “running” of the parameters is as usual determined by the renormalization group 
equation which follows from the behavior of the time ordered product under scaling. 

6 Hamiltonian formalism for quantum field theory, 
and the construction of states 

Up to now we remained in the realm of algebras. There we could study several 
structural properties of the theory. In order to get more detailed predictions of the 
theory one has to evaluate the algebra in specific states. A class of states on the 
local algebras can be obtained in terms of the states of the free theory by embedding 
the interacting theory into the free one, but this is highly ambiguous and gives no 
direct interpretation of the states. Conceptually, one does not need more, since the 
interpretation can be done in terms of the expectation values of observables. In 
practice, however, one would prefer to have states with an a priori interpretation 
as e.g. the vacuum state. The standard way to compute it is the evaluation of 
the product or the time ordered product of interacting fields with an interaction 
Li (go) = f L i (x) go(x)dg,(x) in the vacuum state of the free theory and performing 
the adiabatic limit go —> 1. This limit is well behaved in massive theories, but exists 
also for a suitable sequence (go)n —> 1 hr certain massless theories such as massless 
< p 4 or QED. In the case of time ordered products one just reproduces the standard 
formulas in terms of Feynman graphs; in the case of operator products one has to 
use Steinmann’s sector graphs |42| . The adiabatic limit in this form, however, does 
not always exist, in particular not for states with nonzero temperature. 

A more direct way of constructing states with specific properties could be imag¬ 
ined in a Hamiltonian formalism, as well known from nonrelativistic quantum me¬ 
chanics. The difficulty is that the interaction Hamiltonian for a local QFT is very 
singular so that perturbation theory for selfadjoint operators cannot be used. There 
are two independent reasons for the singular character of perturbations in QFT. 
The first is translation symmetry. In Minkowski space this leads to Haag’s theorem, 
which states that the ground state of the interacting theory cannot be represented 
by a vector in the Fock space of the free theory. If one takes this into account by 
restricting the interaction to a finitely extended spatial region, one can indeed apply 
the perturbation theory of selfadjoint operators in certain superrenormalizable mod¬ 
els in 2 dimensions. One can then construct ground states and consider their limit 
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if the cutoff is removed. In 4 dimensions, however, the local interaction densities are 
too singular, so that also the spatially restricted interaction is not an operator. 

The Hamiltonian formalism relies on a split of spacetime into the product of a 
Cauchy surface and the time axis, and all the observables of the theory are con¬ 
structed in terms of their initial values on this surface, which are supposed to be 
independent of the interaction. But from renormalization theory it is well known 
that in general one has to expect modifications of the canonical structure; moreover, 
even for free fields, the restriction to a Cauchy surface is singular for all nonlinear 
local fields. 

Instead we use the fact that for generic perturbative QFT’s the time-slice axiom 
holds. Moreover, as we saw from the discussion of the proof of this fact, the free and 
the interacting algebra of a time slice can be identified. This suggests to compare 
their time evolutions. Both are automorphism groups acting on the same algebra, 
and they differ by a cocycle. In case of a spatial cutoff of the interaction, the 
cocycle is implemented by a unitary cocycle within the algebra, whose generator is an 
integral over an operator valued function which may be interpreted as a regularized 
interaction Hamiltonian density Kj (x). 

As in section 14.31 we consider the space D n of test functions and the algebra 
generated by S(f) = Ai/’d/z)). We also assume that A 0 = Lj is the 

interaction Lagrangian density. The time slice property proven in section [5] induces 
isomorphisms between the free and the interacting algebras. Let x be a smooth 
function of time t with y(f) = 1 for t > — e and y(f) = 0 for t < — 2e. Then 
supp((t,x) i—^ G(t,x)x(t)) is past compact. We now define a map from 21 G to 21 by 

7x( 5 g(/)) = S Gx {f) , supp / C (—e, e) x K 3 . 

Due to the time slice property this map extends to an isomorphism. Moreover, it only 
slightly changes the kinematical localization at t = 0. Let O r = {(t, x)||i| + |x| < r}. 
Then 

7 x (2l G (0 r )) C 2l(O r+4 e) C 7 x (2l G (0 r+8e )) . 

Let G be constant, let a^’ x = 7 x o a x °7 x -i be the translations of the interacting 
theory mapped to the free theory, and consider the cocycle /3^’ x = a^’ x o a_ x . We 
Gy 

find PtQx) = id anc b f° r / with supp / C O r and small t, 

= Ad S hx (h( X t ~ x))(S(f)) 

where h is time independent, has compact spatial support and h = G on 0 r+ 4 e . 
Proposition 1. The unitaries = Sh x {h(xt — x)) fulfill the cocycle equation 

U%„ = U^ x a t (Ug X ) 

Proof. For sufficiently large u (depending on s, t ) we have 
S h x (KXt - x))ut(S hx (h(x s - X)) = 

s h(x-xu)( h (Xt - X))a t ( s Hx-xu-t)( h (Xs - x)) = 

S h( x -xu)(Kxt - x))Sh( Xt - Xu ){Kxt+ s - Xt)) = 

S(h(x - x„)) _1 S'(/i(xt - Xu))S(h(xt - Xu^Sihixt+a - Xu)) 

= S Hx-xu)( h (Xt+s - x)) = S hx (h(xt+s - X)) ■ 

□ □ 
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We conclude that the unitary cocycle describes the interacting time evo¬ 
lution (with spatial cutoff h) in the interaction picture. Due to the finite speed of 
propagation, it coincides with the full time evolution for small t. 

We now consider a time translation covariant representation (J£,7 t, Uq) and as¬ 
sume that the map 1)9/-) n(S(f)) is strongly continuous. Then the cocycle 
is strongly continuous, and 

t u hx (t) = U? x U 0 (t) (82) 

is a strongly continuous 1-parameter group with selfadjoint generater Hh x which 
describes the dynamics of the interacting system with spatial cutoff. 

In case 7r is irreducible, one may now determine the spectrum of Hh x and inter¬ 
pret it as the energy spectrum of the interacting theory with spatial cutoff (up to 
an additive constant). One may also look for a ground state and consider the limit 
of removal of the cutoff. 

If 7T is a representation induced by a KMS state, and Do is the corresponding 
cyclic vector in the representation space, one knows by Connes’ cocycle theorem that 
there exists a weight whose modular automorphims are the time translations of the 
interacting theory. If Do is in the domain of e~^ Hhx , then this weight is bounded 
and induced by the vector 

V hx = e^ Hh *n 0 . (83) 

If the cocycle is strongly differentiable on a dense domain, the interaction Hamil¬ 
tonian can be defined as the generator of the cocycle. We obtain nu 

H i x = n ^ t u t x = hS{hxT l ^- t S(h X t) =Rv(h X )(Y(hx )), 

where in the last step we have used Bogoliubov’s formula (l55l) for interacting fields. 
In the limit e —> 0, % tends to the ^-function and we obtain the usual interaction 
picture. 

We illustrate the method on the example of an interaction with external sources. 
We start with the CCR algebra 2ls 0 of the free scalar field introduced at the end of 
section 13.51 The formal S-matrix is 

S(f) = S(F f ) = e ^//» e -A</. AD /> , 

where Ff{g>) = f <pfdg. One can verify it by direct computation (using the forumlas 
for time-ordered product given in section 14.21) or, indirectly, by the verification of 
the causal factorization property CO). Namely, we have 

S(f+g)- 1 S(f + g + h)= (84) 

e iFh/h exp % -{(f + g, A D (f + g)) - (/ + g + h, A D (f + g + h)) + (/ + g. Ah)) 

= e^ h) exp A D h) - (f + g, (2A^-A )h)) , 

=a a 

= S(g)- 1 S(g + h)e^’ AA V 

hence if supp /D J_ (supp h) = 0 then by the support property of the advanced 
propagator 

(/, A A h) = 0 
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and the factorization holds. 

We find the interaction Hamiltonian (h time independent) 

Hj X = —tp(hx) — const . 

Due to the smearing in time, this operator remains meaningful also for a pointlike 
source (h ~ <5(x)). 

In general, for the free theory we obtain the usual Fock space Hamiltonian Hq, 
and the Hamiltonian of the interacting theory with spatial cutoff is the sum of the 
free Hamiltonian and the interaction term, 

H = H 0 + j h(x)^Ki(x)d 3 x . (85) 

In this framework, one can now apply the standard perturbative constructions of 
ground states and KMS states. In Bang it was shown that in massive theories 
ground states and KMS states for positive temperatures exist. Some aspects of this 
formalism involving thermal mass were further developped in [14] in conjunction 
with the principle of perturbative agreement |30) . It is hoped that this regularized 
Hamiltonian picture will allow to close the conceptual gap between the standard 
formalism in nonrelativistic quantum mechanics and quantum statistical mechanics 
and the formalism of relativistic QFT. 

7 Conclusions 

We have seen that the concepts of AQFT can be used in renormalized perturbative 
QFT and yield Haag-Kastler nets (in the sense of algebras of formal power series) 
for generic models of QFT. Due to its axiomatic formulation all possible renormal¬ 
ization methods are covered, and one has an a priori characterization of the class 
of renormalized theories associated to a classical Lagrangian, independent of any 
regularization scheme. For practical purposes, it is nevertheless often appropriate 
to introduce a regularization, and in particular analytic regularization schemes such 
as dimensional or analytic renormalization are useful, for computation but also for 
specifying a theory in its class (e.g. by minimal subtraction), see |17j . One may also 
incorporate the ideas of the renormalization flow equation in the sense of Polchinski 
and made rigorous in [32] • This is exposed in pTj . In these notes we restricted our¬ 
selves to scalar field theories. The generalization to other types of field theory have 
been discussed in several papers; fermionic theories can be treated essentially in the 
same way, and gauge theories can be treated after adding auxiliary fields (ghosts 
etc.) and constructing the time ordered products such that BRST symmetry is re¬ 
spected Ea |2 si Eg. Even gravity can be included where however the concept of 
local algebras of observables has to be properly adapted ma- 
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